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ABSTRACT 

We present results from a statistical analysis of 173 bright radio-quiet AGNs selected from the 
Chandra Deep Field-North and Chandra Deep Field-South surveys (hereafter, CDFs) in the redshift 
range of 0.1 < z < 4. We find that the X-ray power-law photon index (P) of radio-quiet AGNs 
is correlated with their 2-10 keV rest-frame X-ray luminosity (Lx) at the > 99.5% confidence level 
in two redshift bins, 0.3 ^ 2; < 0.96, and 1-5 < 2; < 3.3 and is slightly less significant in the 
redshift bin 0.96 ^ 2 < 1.5. The X-ray spectral slope steepens as the X-ray luminosity increases for 
AGNs in the luminosity range 10"*^ to 10''^ erg s~^. Combining our results from the CDFs with those 
from previous studies in the redshift range 1.5 ^ z < 3.3, we find that the P — Lx correlation has 
a null-hypothesis probability of 1.6 xlO^^. We investigate the redshift evolution of the correlation 
between the power-law photon index and the hard X-ray luminosity and find that the slope and offset 
of a linear fit to the correlation change significantly (at the > 99.9% confidence level) between redshift 
bins of 0.3 < z < 0.96 and 1.5 < z < 3.3. We explore physical scenarios explaining the origin 
of this correlation and its possible evolution with redshift in the context of steady corona models 
focusing on its dependency on variations of the properties of the hot corona with redshift. 
Subject headings: cosmology: observations — galaxies: active — galaxies: statistics — X-rays: galaxies 



1. INTRODUCTION 

It is important to extend the study of quasars to high 
redshifts in order to understand their evolution and en- 
vironments. A relevant conclusion from modern studies 
is that the quasar luminosity function evolves positively 
with redshift, having a comovin g space density strongly 
peaked at z pa 2 (e g., ISchmidt lll968t iBovle et al.lll987l : 
I Warren et al.l lT9M). More recent findings suggest that 
the evolution of the space density of AGNs is strongly de- 
pendent on X-ray luminosity (Lx), with the peak space 
density of AGNs mov ing to higher redshifts for more 
lumin ous AGNs (e.g., lUeda et all 120031 : iHasinger et al.l 
[2005h . 

The X-ray band probes the innermost region of the 
central engines of AGNs. The study of AGNs in the 
X-ray band provides important insights about their cen- 
tral engines and the evolution of the AGN luminosity 
function. In most AGNs, the observed X-ray contin- 
uum can be modeled using a power-law of the form 
N{E) = NaiE/Eoy^, where P is the photon index. 
This power-law is attenuated by material in our Galaxy 
as well as material intrinsic to the host galaxy. Sev- 
eral recent studies have centered on estimating the dis- 
tribution of intrinsic column densities (A'h) and the frac- 
tion of AGNs having > 10^^ cm~^. Recent theo- 
retical studies of AGNs (e.g.. iHopkins et ani2005( ) sug- 
gest that the distribution of A^h is luminosity depen- 
dent; this is supported observationally with the detec- 
tion of an anti-correlation between the obscuration frac- 
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tion a nd luminosity (e.g..lSteffen et al.ll2003 : lUeda et al.l 
120031: iLa Franca et all 120051: lAkvlas et all l2006l) . We 

note, however, recent work by iDwellv fc PageT (p006h 
reporting that the obscuration fraction may be indepen- 
dent of luminosity. The dependence of the obscuration 
fraction on redshift is a controversial issue. Some authors 
detect an increase of the obscuration fraction with red- 
shift re.g..lLa Franca et"aI1l2005l : iTreister fc Urrvll2006l : 
iTozzi et al.ll2006D . while others do not find any evidence 
for evolution (e g..lUeda et all l2003l : lAkvlas et all l2006t 
IDwellv fc Pagell2006[ ). 

A recent mini-survey of relatively high-redshift (1.5 < 
z < 4) gravitationally lensed radio-quiet quasa rs (RQQs) 
obser ved with Chandra and XMM-Newton (jDai et al.l 
|2004[ ) indicated a possible correlation between the X-ray 
power-law photon index and X-ray luminosity. This cor- 
relation, characterized by an increase of P with Lx, was 
found for RQQs with 2-10 keV luminosities in the range 
10^'^ to 10^^ erg s~^. Such a correlation is not found in 
nearby z < 0.1 quasars fe.g.. lGeorge et al.ll2000f) . Several 
studies to date of high-redshift quasars do not have large 
enough sample sizes in the 2-10 keV luminosity range 
10^"^ to 10"*^ erg s~^ to place any sig nificant constraints 
on a possible P — Lx c orrelation fe.g.. iReeves fc Turner I 
l200nl: iPage ei~ani2005h . ^ 

One of the concerns with the iDai et al] (|2004l ) analy- 
sis was that the limited number of quasars in the sam- 
ple, combined with the poor signal-to- noise ratio (S/N) 
available for several of the observations and the relatively 
large fraction of BAL quasars, may have led to problem- 
atic systematic effects. Th e number of ava ilable lensed 
radio-quiet quasars used by lDai et al.l ()2004l ) was limited 
to a total of 25 sources, of which the brightest 11 had X- 
ray observations. In order to increase the size of the high- 
redshift radio-quiet quasar sample, we have compiled 
a sample of 173 high-redshift AGNs with moderate-to- 
high S /N spectra available from the Chandra Deep-Field- 
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North and Chandra Deep-Field-South s urveys (CDF-N 
and CDF-S, respectively; i ointly CDFs; iGiacconi et al.1 
[2OOI [Alexander et al.l[200l . 

The main scientific goal of this work is to constrain 
better the F — Lx correlation found by Dai et al. (2004). 
The significant increase in sample size allows us to place 
tighter constraints on the significance of the correlation. 
We also test the correlation in narrower redshift bands 
which will allow us to determine the epoch by which pos- 
sible changes in the average emission properties of AGNs 
occurred. Currently, the two deepest X-ray surveys are 
the CDF-N and CDF-S with «2 Ms and wl Ms expo- 
sures, respectively. Both surveys cover wSOO arcmin^ ar- 
eas and target different regions of the sky characterized 
by low Galactic column densities and an absence of bright 
stars ( Giacconi et al. 2002; Alexandcr_et_al- 2003). The 
CDFs pointings have sufficient sensitivity to detect the 
X-ray emission from AGNs with moderate luminosities 
(Lx « 10"^ - 10'*'' erg s"*) out to z « 2 - 6. 

Radio-quiet AGNs (RQ AGNs) correspond to the ma- 
jority of active galaxies (~ 90%) that contain a central 
active nucleus and show several differences in their spec- 
tral properties compared to radio-loud AGNs. Radio- 
loud AGNs have powerful sub-parsec jet-linked X-ray 
synchrotron self-Compton (SSC) emission, which intro- 
duces an additional component to their spectra. As a 
consequence RQ AGNs are observed to have, on av- 
erage , steeper X-ray pow er-laws than radio-loud AGNs 
(e.g.. iReeves et allll997[ ). We therefore have chosen to 
exclude radio-loud AGNs from this study. Throughout 
this paper we adopt a flat A-dominated universe with 
Ho=70 km s-iMpc^\ I^a = 0.7, and f^M = 0.3. The 
Chandra data were reduced using the CIAO version 3.3 
software tools provided by the Chandra X-ray Center 
(CXC), and the spectral analysis was performed using 
XSPEC version 12. 

2. SAMPLE SELECTION 

Our sources were selected from the CDFs, cur- 
rently the two deepest X-ray surveys. The on- 
axis sensitivity limits for the CDF-N are « 2.5 x 
10-1^ erg cm-2 s'^ (0.5-2.0 keV) and « 1.4 x lO'^^ 
erg cm~^ (2-8 keV). These limits are around two 
times more sensitiv e than those for the CDF-S 
(jGiacconi et al.l l2002t [Alexander et al.l l2003f ). The 
CDFs are 50-250 times more sensitive than previ- 
ous X-ray surveys, detecting « 900 point sources, of 
which « 600 are AGNs and gala xies with measured 
redshifts (IGiacconi et all [20 02; Al exander et al.l 120031 ; 
iBarger et al.ll2003l ; IZhengeral.ii2004i r 

Spectroscopic and photometr ic redshifts were 

gathered fro m the literature (ICroom et al.l l2001l; 

Barger et ahl 120031; ISteidel et al l 120031; ICowie et al l 



2004; Mobasher et al. 2004 
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Wirth et al 

2004 lAle'xand cr et al 



Wolf et al.^ 
2005 



" iSzokolv et al 
'2004!; 



2004 



_ _ ■ ■ _ _ _ , Colbe rt et al . . 

Le Fevre et al.ll2005l; IVanzella et al.ll2006[ ) and vetted to 



Zheng et al.l 
[200a 



remove redshifts which did not appear to belong to the 
most-likely optical counterpart to each X-ray source. 
The latter was assessed by comparing the optical and 
X-ray images, which were aligned to betwee n 0'.^12-0"25; 
this n otably affected the redshifts from IZheng et al.l 
()2004l ). where 47 («14%) of the redshifts were rejected 
for being associated with an unlikely optical counterpart. 
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Fig. 1. — (upper panel) Number of sources with more than 
S photons (0.5-8 keV observed-frame) versus S. The thick hue 
corresponds to all CDFs sources with measured redshifts of ^ > 0.1. 
The dotted line shows sources of the CDF-N survey with z '> 0.1, 
and the dashed line shows sources of the CDF-S survey with z > 
0.1. The vertical dotted line corresponds to 5=170. Note that the 
sample used to generate this figure contains AGNs (both RQ and 
radio-loud AGNs), normal galaxies and starburst galaxies, (lower 
panel) Number of radio-quiet AGNs with z > 0.1 vs. the number 
of photons {S; 0.5—8 keV) in their spectra. The solid line represents 
sources with fits performed in the 0.5—8 keV observed- frame band 
and the dotted line represents sources with fits performed in the 
2-10 keV rest-frame band. 



For the ^^40 faint sources still lacking redshift estimates, 
we used the BPZ code (B ayesian photometric redshift 
estimation; iBenitez I [2OOOI) and available photometry 



(lArnouts et al.ir2001l ; ]Barger et al.ll2003l;lGiavansco et all 
|2004|) to estimate crude redshifts. 

The selection criteria for our sample of RQ AGNs are 
(1) that the sources are radio-quiet (see below), (2) that 
the redshifts of the sources are greater than 0.1, and (3) 
that the total number of photons in the full band (0.5-8 
keV) is greater than ^ 170 counts {S > 170) result- 
ing in moderate-to-high S/N spectra. The selection of 
a cut-off at ~ 170 counts allows an accurate estimate 
of the photon index, wh ich is not possible for fainter 
sources (jTozzi et al.ll2006[) . Based on the condition that 
S > 170, the on-axis flux limits of our sample in the full 
band (0.5-8 keV) in the CDF-N and CDF-S surveys are 
« lxl0~*^ erg cm^^ s^^ and ~ 2xl0~^^ erg cm^^ s^* , 
respectively. 

In Figure [1] (upper panel), we show the cumulative dis- 
tribution for number of X-ray sources having more than 
S counts (0.5-8 keV) for the CDF-N and CDF-S. The 
CDFs contain 205 sources with more than 170 counts 
at z > 0.1. Most of these sources are AGNs; however, 
in the low-redshift regime of our sample 0.1 < z < 1.0 
we expect onl y a small fraction of starb urst and "nor- 
mal" galaxies (jBrandt fc Hasingerll2005f l. Following the 
classification scheme discussed in §4.1.1 of ^ Bauer et al.l 
(|2b04a1 . we found two starburst galaxies and one "nor- 
mal" galaxy, which we remove leaving 202 AGNs in our 
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sample. 

Radio-loud AGNs were classified based on a radio- 
loudness parameter R > 10 {R = /sghz/Zb)- To find 
these sources, we matched the X-ray positions with ra- 
dio sources using a matching radius of 2 arcsec. The 
flux-density at 5 GHz was obtained from the flux-density 
at 1.4 GHz assuming a power law radio spectrum {f^, oc 
i/^"''), where = 0.8 is a characteristic radio spectral 
index of synchrotron radiation^. 

The flux in the B filter was obtained from Bar ge r et all 
|2003) for the CDF-N sources and from public-domain 
tables of the GOODS and COMBO-17 surveys for the 
CDF-S sourc es. When searching the radio catalogs^ 
provided by iRichards I (|2000[ ) for the CDF-N and 
lAfonso et all (12006 1) for the CDF-S, we find that 29 
(~14%) out of the 202 X-ray detected AGNs were radio- 
loud. This leaves 173 RQ AGNs which we use for our 
analysis out of which 111 have spectroscopic redshifts. 

3. SPECTRAL EXTRACTION 

The X-ray spectra of the sources of the CDFs analyzed 
in our study were extracted using the software routine 
ACIS_EXTRACT v3.94 (hereafter AE; Townsley et al. 2003; 
Broos et al. 2005), included in the Tools for ACIS Real- 
time Analysis (TARA; Oct 20, 2005) software package. 
^ AE is ideal for extracting and analyzing the spectra of 
large numbers of point and diffuse sources observed with 
ACIS over multiple epochs. AE calls procedures from 
both CIAO (v3.3) and HEASOFT (v6.0.4) and uses cal- 
ibration files that are part of the CALDB v3.2.1 product 
provided by the Chandra X-ray Center. 

The w2 Ms CDF-N (f«l Ms CDF-S) observations are 
comprised of 20 (10) event files. The event files were 
corrected for charge transfer inefficiency, bad columns, 
bad pixels, and cosmic ray afterglows. The event files 
were also filtered for time intervals of acceptable as- 
pect solution and background levels. A detailed descrip- 
tion of the data reduc tion procedures are presented in 
[Alexander et al.l (|2003f ). Background event files and ex- 
posure maps were created by excluding circular regions 
centered on the detected sources with radii that are a 
factor of 1.1 times larger than the 99% encircled energy 
radii of the point spread functions at ^ 1.49 keV. Source 
extraction regions were constructed to contain 90% of the 
PSF encircled energy derived from the CXC 1.4967 keV 
PSF libraries. There were two exceptions to this proce- 
dur e. First, for sources with greater than 1000 counts in 
the [Alexander et al.l (|2003l ) catalog we used extraction 
regions that contained 99% of the PSF encircled energy. 
Second, for sources with 90% encircled energy extrac- 
tion regions that overlapped we reduced the extraction 

^ In A GNs values of ar could be flatter than the adopted ar = 
0.8 fe.g.. lRlchards et al . 1998; Muxlow ct al. 2005), with measured 
standard deviations ~1 (e.g., Wadadekar 2004). We investigated 
how a flatter ar may affect our results and find that choosing a 
value of ar = 0.6, for example, to estimate R will not change our 
sample of RQ AGNs, whereas, a value of ar = 0.4 will result in 
the exclusion of only two sources from our sample (1% of the entire 
sample) in order to satisfy R < 10. We conclude that our sample 
selection and results of our statistical analysis are not significantly 
affected by values of the radio spectral index as low as ar = 0. 

6 The radio surveys of Richards (2003) and lAfonso et al.l pOOeD 
cover the entire CDF-N and CDF-S regions respec tively. More 
details of the CDF-S radio observations are found in lNorris et al.l 
igOOi). 

TARA is available at |http : //www ■ astro ■ psu ■ edu/xray/docs/TARA/ 1 



regions to avoid overlap. Local background extraction 
regions were chosen as annuli centered on the source po- 
sitions with inner radii equal to that of the source ex- 
traction regions and with outer radii selected such that 
the background region contained at least 100 background 
counts and had an area at least fours times that of the 
source region. 

We note that in the current analysis we made no at- 
tempt to correct for possible spectral variability over the 
few year period of the observations of the CDFs. Spectra 
obtained are therefore time-averaged over the period of 
the observations. 

4. SPECTRAL ANALYSIS 

Two energy bands were used to fit the Chandra spec- 
tra: the 0.5-8 keV observed-frame and the 2-10 keV 
rest- frame. To obtain the maximum S/N we utilized the 
observed-frame energy range of 0.5-8 keV. The lower en- 
ergy bound was chosen because the Chandra effective 
area is not well calibrated below 0.5 keV, and the up- 
per energy bound was chosen because the S/N decreases 
greatly above this energy for most of the sources in the 
sample. One advantage of using the same observed-frame 
energy range for every object is that the same systematic 
instrumental uncertainties apply to every fit. Since most 
of the detected spectrum is used in the analysis, the S/N 
is higher than for cases where restricted energy ranges 
were used. 

To test how the F — Lx correlation might be affected by 
absorption and possible contamination from other emis- 
sion processes, we also fitted the spectra in the rest-frame 
energy range of 2-10 keV. This range was selected to 
avoid possible contamination from soft-excess emission 
that is often detected in AGNs below rest-frame energies 
of ^ 1 keV. The selection of the 2-10 keV rest-frame 
band also aids in reducing the effects of X-ray absorp- 
tion. For example, assuming a source with a power-law 
spectrum of r=1.7, z=l, A'n^lO^^cm"^, and solar abun- 
dances the fraction of absorbed photons is 30% in the 
0.5-8 keV observed-frame band and 9% in the 2-10 keV 
rest-frameband. The 2-10 keV rest-frame also mini- 
mizes possible contamination from Compton-reflection 
emission from circumnuclear material that is thought to 
peak at a rest-frame energy of about 20 keV. In gen- 
eral, 2-10 keV rest-frame spectra have fewer counts than 
0.5-8 keV observed- frame spectra. For fits performed in 
the 2-10 keV rest-frame band, we selected sources with 
more than 170 counts in this band, leaving a sub-sample 
of 144 RQ AGNs. 

The total number of photon counts per source {S) 
with energies in the 0.5-8 keV observed-frame band lies 
in the range 170-13000. In Figure [1] (lower panel) we 
present the number of z > 0.1 radio-quiet AGNs in 
our sample versus the number of photons with energies 
in the 0.5-8 keV observed-frame band. The solid line 
applies to sources with spectral fits performed in the 
0.5-8 keV observed- frame, and the dotted line applies 
to sources with spectral fits performed in the 2-10 keV 
rest-frame. The mean logarithm of S for sources with 
spectral fits performed in the 0.5-8 keV observed- frame 
is (log 5')=2.74 with a standard deviation of cr ~ 0.42. 
The mean logarithm of S for sources with spectral fits 
performed in the 2-10 keV rest-frame is (log S')=2.83 
with a standard deviation of cr 2± 0.41. Based on the 
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TABLE 1 

Models used in fitting the spectra of the RQ AGNs of 
our sample.^ 



50 f 



Model b 


No sources 


% of the whole sample 


PL 


77 


44.5 


APL 


76 


43.9 


PAPL 


9 


5.2 


lAPL 


4 


2.3 


PL+EL 


4 


2.3 


APL+EL 


3 


1.7 



^ The selection criteria for the sample of RQ AGNs of our 
present study were that the redshifts of the sources were 
greater than 0.1 and the total number of photons in the full 
band (0.5-8 keV) was greater than ~ 170. ^ PL=power-law 
(XSPEC model wabs(pow)); APL=absorbed power-law (XSPEC 
model wabs*zwabs(pow)); PAPL=partially absorbed power-law 
(XSPEC model wabs*zpcfabs(pow)); lAPL =ionized absorbed 
power-law (XSPEC model wabs*absori(pow)); PL-f-EL^power- 
law + emission-line (XSPEC model wabs(pow+zgauss)); 
APL-t-EL=absorbed power-law -t- emission-line (XSPEC model 
wabs*zwabs* (pow-t-zgauss) ) . 

fact that our sample contains so urces with re latively low 
counts, we used the C-statistic (|Cash Ill979t ) to fit spec- 
tra as adopted in a similar study presented in lTozzi et al.l 
([2006'). In this study, the authors concluded that the 
C-statistic is more accurate than the x^-statistic in esti- 
mating the spectral parameters of AGNs with low-count 
sp ectra (~ 100 counts); si milar arguments are presented 
in iNousek &: Shuel (|1989| ). We also performed spectral 
fits in the 0.5-8 keV observed-frame band using the 
statistic, with a grouping of 10 counts per bin. The sole 
purpose of using the x^-statistic was to apply the F-tcst 
to assess the use of more complex spectral models. 
For the CDF-S and CDF-N sources of our sample, we 
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assumed Galactic column densities of 8^8 

(jStark et al.l(T993 ) and 1.3x10^0 cm'^ (iLockman Il200l . 
respectively. The spectral analysis was performed us- 
ing XSPEC version 12. The default spectral model 
used is a power law (PL; POW) with Galactic ab- 
sorption (WABS). Additional model components were 
added to the default model in cases where the F-test 
showed an improvement in the fit at the 95% confi- 
dence level (0.5-8 keV observed-frame) when these addi- 
tional components were used. We refer to models com- 
prised of the default model plus additional model com- 
ponents as alternative models. Alternative models in- 
cluded an absorbed-power-law model (APL) at the red- 
shift of the source (WABS ZWABS POW), an ionized- 
absorbed-power-law model (lAPL) (WABS ABSORI 
POW), a partial-absorbed-power-law model (PAPL) 
(WABS ZPCFABS POW) and/or models that included 
an iron line (PL-hEL;APL-hEL)(WABS ZGAUSS POW; 
WABS ZWABS ZGAUSS POW). We also considered 
models that contained a Compton-reflection component 
(PEXRAV) , but we did not find any improvement in the 
fits using such models. Our finding, that none of the 
sources in our sample require a Compton -reflection com- 
ponent, is in agreement with lTozzi et al.l 12006,) who find 
that only 14 out of 321 CDF-S sources require Compton- 
reflection components. We note that none of these 14 
sources are part of our sample mostly because they con- 
tain less than 170 counts. Even though we do not detect 
a significant refiection-componcnt in our spectral fits, an 
unaccounted reflection-component could still be affecting 
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Fig. 2. — Number of sources vs. luminosity (upper panel) 
and redshift (lower panel). The solid line represents sources with 
fits performed in the 0.5-8 keV observed-frame band and the dot- 
ted line represents sources with fits performed in the 2-10 keV 
rest-frame band. 

the estimation of F (see §5.2.6 for details). In Table [1] 
we list the number of sources from our sample (the en- 
tire sample contains 173 RQ AGNs) fit with a particular 
spectral model. From Table [2 we notice that based on 
the F-test there are 92 sources with detected absorption 
(~53 % of the whole sample) and 7 cases with detected 
iron lines (~4 % of the whole sample). 

The spectral-fitting results are presented in Table [2] 
for fits performed in the 0.5-8 keV observed-frame band 
(173 RQ AGNs), and the fits performed in the 2-10 keV 
rest-frame band (144 RQ AGNs). In Table H we pro- 
vide the photon index F (errors at the 90% confidence 
level), the intrinsic column density iVn (errors at the 
90% confidence level) in units of 10^^ cm~^, and the log- 
arithm of the hard X-ray luminosity in the rest-frame 
2-10 keV band in units of erg s~^ (hereafter referred to 
as L2_io). Table[2]also includes the X-ray identification 
of the sources based on their RA and DEC positions, the 
photon counts in the fitted range, the number of degrees 
of freedom, and the values of the C-statistic. The last 
two quantities provide an estimate of the quality of the 
fits. 

In Figure [21 we show the distributions of 2-10 keV 
luminosities (upper panel) and redshifts (lower panel) 
of the sources in our sample with fits performed in 
the 0.5-8 keV observed-frame band (solid fine; 173 RQ 
AGNs) and the 2-10 keV rest-frame band (dotted line; 
144 RQ AGNs). The luminosities of the sources in our 
sample cover the range 3 x 10^^ — 6 x 10'*'' erg s~^ where 
the lower limit is mostly determined by the sensitivity 
limit of the CDF-N survey, while the upper limit is a 
statistical consequence of the fact that luminous AGNs 



(^2-1 



> 



lO*'^ erg s^i) 



are less numerous than lower 



luminosity AGNs (see e.g., iBrandt &: Hasinger 1 120051. 
The mean redshift and mean logarithmic X-ray luminos- 
ity of the sources with fits performed in the 0.5-8 keV 
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TABLE 2 

Properties of our, sample of RQ AGNs selected from the Chandra Deep Field Surveys. 



Xid z^^ Counts V Nn"^ log L2-10 C-stat dof type<= model ' 



RESULTS BASED ON FITS PERFORMED IN THE 0.5-8 KEV OBSERVED-FRAME. 
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RESULTS BASED ON FITS PERFORMED IN THE 2-10 KEV REST-FRAME. 
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Note. — Table 2 is presented in its entirety in the electronic edition of the Astronomical Journal. A portion is shown here for guidance 
regarding its form and content. 

^ Xid with RA+DEC coordinates.^ Spectroscopic (sp) and photometric (ph) redshifts gathered from the literature (see §2).*^ Back- 
ground subtracted source counts in the 0.5—8 keV observed- frame band for fits performed in the 0.5-8 keV observed-frame band 
and counts in the 2—10 keV rest-frame band for fits performed in the 2—10 keV rest-frameband."^ In units of lO'^^ cm~^.^ Based 
on I Bauer et al] ll2004aD . source classifications from §4.1.1 (http://www.astro.psu.edu/~nicI/hdf/hdf-chandra.html) ^ PL=power-law 
(XSPEC model wabs(pow)); APL=absorbcd-powcr-law (XSPEC model wabs*zwabs(pow)); PAPL=partially-absorbed-power-law (XSPEC 
model wabs*zpcfabs(pow)); lAPL =ionized-absorbed-power-law (XSPEC model wabs*absori(pow)); PL+EL^power-law + emission-line 
(XSPEC model wabs(pow+zgauss)); APL+EL^absorbed-power-law + emission-line (XSPEC model wabs*zwabs*(pow+zgauss)).S The 
estimated values of F, A^h, and X-ray luminosities presented in this table are based on spectral fits performed in the 0.5—8 keV 
observed- frame. 



observed-frame band are (z) ~ 1.41 and (log L2-10) — 
43.6, respectively. The mean redshift and mean loga- 
rithmic X-ray luminosity of the sources with fits per- 
formed in the 2-10 keV rest-frame band are (z) ~ 1.38 
and (log L2_io) — 43.6, respectively. 

In Figure [3] (upper panel), we show the distributions 
of the photon indices of the sources with fits performed 
in the 0.5-8 keV observed-frame band (solid line) and 
the 2-10 keV rest- frame band (dotted line). We find 
the mean photon indices and their standard deviations 
for fits performed in the 0.5-8 keV observed- frame and 
2-10 keV rest-frame bands to be (F) ~ 1.60 ± 0.27 and 
(r) ~ 1.70 ± 0.29, respectively. In Figure [3] (lower 
panel), we show the distributions of the intrinsic col- 
umn densities of sources with significant absorption (only 
sources where the i^-test indicated significant intrinsic 
absorption at the >95% confidence level are included 
in the distributions); sources with fits performed in the 
0.5-8 keV observed- frame band are indicated with the 



solid line (92/173), and sources with fits performed in the 
2-10 keV rest-frame band are indicated with the dotted 
hue (82/144). In the two fitted energy ranges, the peak 
of intrinsic column density distribution is log iVn ^ 22.6, 
and there is a fraction of ^^40% sources from the to- 
tal sample having log A^h > 22. We note that there 
are likely systematic errors on these column density es- 
timates due to unmodeled absorption complexity. These 
parameter valu es and distributions are i n agreement with 
those found in iTozzi et al.l (|2006( l and iDwellv fc Pajel 
('2006^. 

In Figure m we present a diagram comparing estimates 
of F obtained from fits performed in the 2-10 keV rest- 
frame (Frost) and fits in the 0.5-8 keV observed- frame 
(robs)- The size of each symbol in Figure |4] increases with 
redshift. Deviations from the straight line (Fobs=rrest) 
are most likely statistical in nature, however, a few may 
be associated with the effects of intrinsic absorption, soft 
excesses, non-detected spectral lines, and Compton re- 
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Fig. 3. — Number of sources vs. photon index (upper panel) 
and column density (lower panel). In the lower panel we show 
only the sources with a significant detection of absorption in their 
spectra. The solid line represents sources with fits performed in 
the 0.5-8 keV observed-frame band and the dotted line represents 
sources with fits performed in the 2-10 keV rest-frame band. 

flection. 

In general, the agreement between r^est a-nd Fobs is 
good; this is first quantified by a high Pearson linear cor- 
relation coefhcient (^--^0.73) and a very low null hypoth- 
esis probability (~ 4.8 x 10~^^). Secondly, this agree- 
ment is quantified by testing whether the linear relation 
between T^est and Fobs is consistent with Fpest = Fobs- 
To verify the later we performed a fit to the data 
assuming Fj-ost — Q^Fobsi where a was a free parame- 
ter.^ We considered the errors in both variables Fj-cst 
and Fobs when performing the least-squares fit. We 
obtained a = 0.996 ± 0.08 (error at the 68% confi- 
dence level) with ~ 99.3 for 143 degrees of freedom 
(dof). As a basic check for the luminosity dependence 
of the linear relation between F,.ost and Fobs we per- 
formed fits of the model F^ost = «Fobs to sources with 
log L2-10 ^ 43.6 and sources with log L2-10 ^ 43.6. 
We obtained a = 0.995 ± 0.011 (x^ = 43.8; dof=70) for 
sources with log L2-10 < 43.6 and a = 0.997 ± 0.010 
(x^ = 55.6; dof=72) for sources with log L2-10 > 43.6. 

5. RESULTS AND DISCUSSION 
5.1. Selection Effects. 

We begin this section by describing several selection 
effects that are intrinsic to our sample. Having an ex- 
posure of w 2 times greater than the CDF-S survey, the 
CDF-N survey drives the sensitivity limits of our sources, 
so the following discussion will be focused on this survey. 

Figure [5] shows the estimated best-fit column density 
versus the 2-10 keV luminosity (upper panel), and the 

* A fit using the relation y=Qx to model some bivariate sam- 
ple {xi, yi) with errors in both variables {(Jxi, o-yi) is obtained by 
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Fig. 5. — Estimated best-fit column densities versus 2—10 keV 
luminosities (upper panel), and 2-10 keV luminosities versus red- 
shifts (lower panel) of the 2 > 0.1 RQ AGNs. Filled circles repre- 
sent sources with fits performed in the 0.5-8 keV observed-frame 
band and open squares represent sources with fits performed in 
the 2-10 keV rest-frame band. In the upper panel, the two lines 
indicate the maximum column density that can be detected for 
a source with 0.1 < 2 < 4.0, and a total of 170 counts in the 
0.5—8 keV observed- frame (dashed) and 2-10 keV rest-frame (dot- 
ted). In the lower panel the lines indicate the minimum luminosity 
required for the detection of a source as a function of redshift. A 
total of 170 counts in the 0.5-8 keV observed-frame band (dashed) 
or a total of 170 counts in the 2-10 keV rest-frame band (dotted) 
is assumed. For the limits shown in this figure it is assumed that 
the source is detected at the ACIS-I aim-point with an exposure 
time of Ri2 Ms. 
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Fig. 6. — Spearman correlation coefficients of the F — L2-10 
relation as a function of the mean redshift within each sub-sample. 
Each sub-sample contains 38 RQ AGNs. The solid line corresponds 
to the fits performed in the 0.5-8 keV observed-frame band (Table 
[J}. The dashed-line corresponds to the fits performed in the 2— 
10 keV rest-frame band (Table |2]l. The dotted lines correspond to 
three different levels of significance (68%, 90% and 99%), assuming 
38 independent measurements of F vs. L2— lo- 

2-10 keV luminosity versus redshift (lower panel), for 
our z > 0.1 RQ AGNs. Spectral fits performed in the 
0.5-8 keV observed-frame band and 2-10 keV rest-frame 
band are shown with filled circles and open squares re- 
spectively. In the upper panel of Figure[5l the dashed line 
shows the maximum column density that can be found 
for a source with 0.1 < z < 4.0 assuming the source 
is detected at the ACIS-I aim- point with an exposure 
time of w2 Ms ()Alexander et al.lf200a ) and a total of 170 
counts in the 0.5-8 keV observed-frame band (dashed 
line) and 2-10 keV rest-frame band (dotted line). Each 
point on these curves is obtained by fixing A'h and find- 
ing the minimum luminosity that can be obtained with 
0.1 < z < 4.0 assuming a source with F — 1.6, Galac- 
tic column density of 1.3xl0^°cm~^ and 170 counts in 
each fitted energy range. For low-luminosity sources 
(L2-10 ^ 10''^ erg s^^), the threshold column density 
is «10^^ cm~^, which increases by a factor of '^lO for 
higher luminosity sources (lO^'' — 10**^ erg s^^). The 
maximum column density observed at a specific lumi- 
nosity is set by AGNs with z = 0.1 in most of the ob- 
served luminosity range; however, for luminous sources 
(L2-10 > 10'^^ erg s-i) higher redshift AGNs (z ~ 1-4) 
establish the limit in Nh because these are less affected 
by absorption. 

In the lower panel of Figure [5l the curves indicate 
the minimum luminosity required for the detection of 
a source as a function of redshift. We have assumed 
a source free of intrinsic absorption, positioned at the 
ACIS-I aim-point with an exposure time of «2 Ms, 
F = 1.6, galactic column density of LSxlO^^cm"^ and 
with 170 photon counts in the 0.5-8 keV observed-frame 
(dashed line) or 170 counts in the 2-10 keV rest-frame 
(dotted line). The threshold luminosity is «10''^ erg s^^ 
for z « 0.5 and wSxlO'*^ erg s~^ for z ~ 2.5. The dashed 
curve in Figure [5] is obtained by assuming no intrinsic 
absorption; how ever, the presence of A^h, which might 
be evolving fe.g..'La Fr anca et al.ll2005l : iTreister fc Urrv"! 
[2006; Tozzi ct al. 2006), could be increasing the observed 
threshold luminosity. 
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Fig. 7. — Spearman correlation coefficients of the F — L2-10 
relation as a function of redshift for the RQ AGNs within each in- 
dependent redshift bin. In the upper panel each redshift bin con- 
tains ~55 RQ AGNs and the fits were performed in the 0.5-8 keV 
observed-frame band (Table [2J ■ In the lower panel each redshift 
bin contains ~45 RQ AGNs, and the fits were performed in the 
2-10 koV rest-frame band (Table [2]l. The dotted lines correspond 
to three different levels of significance (68%, 90% and 99%); these 
are obtained assuming 55 sources in the upper panel and 45 sources 
in the lower panel. 

5.2. Luminosity and Photon index 

One of the goals of this work is to examine a possible 
correlation between Lx and F in a sample o f RQ A GNs 
which was previously reported bv iDai et al.l (120041 ). To 
improve on the Dai et al. analysis we significantly in- 
creased the sample size using the CDFs, considered a 
larger redshift range, and used X-ray spectra that con- 
tained more than 170 counts in the full band (0.5-8 keV). 
The results of this analysis is shown in the following sec- 
tions. For the following analysis we use the X-ray lumi- 
nosity in the 2-10 keV rest-frame band (L2-io)- 

5.2.1. Possible evolution of ttie strength and significance of 
the F — Lx correlation 

As a first approach, we searched for a F - Lx corre- 
lation as a function of redshift selecting sub-samples or- 
dered in redshift. We used sub-samples containing 38 
sources for fits performed in the 0.5-8 keV observed- 
frame band and the 2-10 keV rest-frame band (Tabled). 
We calculated the mean redshift of each sub-sample and 
computed the Spearman rank correlation coefficient and 
significance of the correlation between F and £2-10. This 
process was repeated by shifting the sampling window 
across the entire observed redshift range. Figure [H] shows 
the values of the significance of the correlations and the 
Spearman correlation coefficients of the F — L2-10 rela- 
tion as a function of the mean redshift of the sources 
within each sub-sample, using the best-fit parameters 
from Table [2] The solid line in Figure [6] corresponds 
to sources fitted in the 0.5-8 keV observed-frame band 
and the dashed line corresponds to those fitted in the 2- 
10 keV rest-frame band. We notice that the correlation 
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TABLE 3 
Correlation table of F versus Lx- 
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^ Number of RQ AGNs in each redshift 
correlation.'^ Calculated from T versus log Lx- 
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Fig. 8. — F versus 2-10 keV luminosity of radio-quiet AGNs in the 
redshift range of 0.3 ^ 2 < 3.3. In the upper panel we show sources 
with fits performed in the 0.5—8 keV observed-frame band. In the 
lower panel we show sources with fits performed in the 2-10 keV 
rest-frame band. The symbol size increases with redshift. Filled 
squares are sources with 0.3 ^ 2 < 0.96, open circles are sources 
with 0.96 ^ 2 < 1.5, and filled triangles are sources with 1.5 < 
2 < 3.3. The dotted line indicates the least-squares fit to sources 
having 0.3 ^ 2 < 0.96. The dashed line shows the least-squares fit 
to sources having 1.5 ^ 2 < 3.3. 

has two significant peaks (~99%) in both energy bands 
fitted, one with a mean redshift of ~ 0.7, and the other 
with a mean redshift of ~ 2.2. 

As a second approach, we selected three independent 
redshift bins covering the redshift range 0.3 < z < 3.3. 
The high redshift bin (1.5 < z < 3.3) w a s cho sen to 
match the redshift range where iDai et al] (|2004f ) found 
the r — Lx correlation while the other two redshift bins 



(0-3 <z< 0.96 and 0.96 <z< 1.5) were selected to ob- 
tain independent redshift bins with comparable numbers 
of sources within them. Each redshift bin contained ~ 55 
sources in the 0.5-8 keV observed-frame band and ^ 45 
sources in the 2-10 keV rest-frame band. In Figure [7] 
and Table [3] we show the Spearman rank correlation co- 
efficient and the significance of the Spearman correlation 
coefficient in each bin. The upper panel of Figure [7] corre- 
sponds to fits performed in the 0.5-8 keV observed-frame 
and the lower panel to fits performed in the 2-10 keV 
rest-frame. The height of each bar is the significance of 
the r — ^2-10 Spearman correlation. The correlation for 
fits performed in the 0.5-8 keV observed-frame is signifi- 
cant for the three redshift bins; however, we find a slight 
decrease in the strength and significance in the second 
redshift bin for fits performed in the 2-10 keV rest-frame. 
The significance of the correlation in the first and third 
redshift bins is >99.5% for both fitting ranges. A sig- 
nificant expansion of our sample made by incorporating 
additional deep AGN surveys will be required to confirm 
the possible decrease of the strength of the correlation in 
the second redshift bin. 

In FigurelH we plot F vs. ^2-10 for sources in each red- 
shift bin of Figure [3 for fits performed in the 0.5-8 keV 
observed-frame (upper panel) and the 2-10 keV rest- 
frame (lower panel). Sources in our sample with redshifts 
in the range 0.3 ^ z ^ 0.96 have a lower mean luminosity 
of (log L2-10) 43.1 (cTiogia-io ~ 0-5) than sources with 
redshifts in the range 1.5 < z < 3.3 which have a mean 
luminosity of (log L2-10) ^ 44.1 (fiiogLa-io ~ 0-4)- The 
luminosity distributions for sources in the redshift bins 
0-3 ^ z < 0.96 and 1.5 < z < 3.3 are shown in Figure O 
We notice that the peak of the distribution of the sources 
in the high-redshift bin is significantly higher in luminos- 
ity than the peak of the distribution of the sources in 
the low-redshift bin. This shift in luminosity distribu- 
tions is mainly a selection effect (see Figure [5]) combined 
with the fact th at luminous sources are more numerou s 
at high redshift (|Ueda et al.ll2003l : iHasinger et al.ll2005[) . 
The Spearman correlation index of the F — L2-10 data 



9 



20 \ 



15 ■ 



0.3 S z S 0.96 
1.5 S z S 3.3 



E 10 - 



5 - 



OLE 



41.5 42.0 42.5 45.0 45.5 44.0 44.5 45.0 
log L2_,;, 

Fig. 9. — 2—10 keV rest-framo luminosity (L2— lo) distributions 
for radio-quiet AGNs with 0.3 ^ 2 ^ 1-5 (thick hne) and 1.5 < 
z < 3.3 (dashed Une). The fits are performed in the 0.5-8 keV 
observed-frame band. 

for the whole sample (173 RQ AGN) is -0.24 (99.8% 
significance) and ^^0.16 (94.1% of significance) for the 
0.5-8 keV observed-frame and the 2-10 keV rest-frame 
spectral fits. These correlation coefficients are signifi- 
cantly lower than those found in the 0.3 z < 0.96 and 
1-5 z < 3.3 redshift bins (see Table[3]for more details). 

In Table [H we show the results of a test of the F — Lx 
correlation using only sources with spectroscopic red- 
shifts. We find that the F — Lx correlation of the sub- 
sample of sources with spectroscopic redshifts is signif- 
icant in the three redshift bins; however, as indicated 
in Table m this sub-sample includes a larger fraction of 
type 1 AGNs and contains more sources with log A'h^ 
22 than that of the whole sample. In addition, the size 
of this sub-sample is significantly smaller than the whole 
sample. We caution that the strengths of the correlations 
provided by the non-parametric tests used in our analysis 
of sub-samples containing a small number of sources N 
with the present uncertainties in the photon indices may 
be inaccurate since the variance of the Spearman corre- 
lation coefficient is 0^=^^^. Photometric redshifts are 
subject to larger errors than spectroscopic ones and for 
sources with z > 1 the error is approxim ately given by 
Az/(1 + z) = 0.05 fe.g. lCohen et al.ll2000[ ). In our anal- 
ysis, the uncertainties in the redshifts will mainly affect 
the estimation of the X-ray luminosities. For example, a 
source at 2: ~ 2 will have an uncertainty in the estimated 
luminosity of ALx/Lx ^ 0.3. This level of uncertainty 
will not significantly affect our results since our study in- 
volves estimating changes in the photon index over two 
orders of magnitude in X-ray luminosity. 

In the following three sections, we focus on sources 
in the first bin (0.3 ^ z < 0.96) and the third bin 
(1-5 ^ z < 3.3) and test the sensitivity of the F — Lx 
correlation to the possible presence of intrinsic absorp- 
tion and Compton reflection in the spectra of the sources. 

5.2.2. Possible evolution of the slope and offset of the 
r — Lx correlation 

In Figure [TUl we show F versus L2_io for sources in the 
ranges of 0.3 < z < 0.96 (upper panel), and 1.5 < z < 3.3 
(lower panel). The values of the X-ray luminosities and 
spectral indices shown in Figure [TUl were obtained by fit- 
ting the spectra in the observed-frame energy range of 
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Fig. 10. — r versus 2-10 keV luminosity of radio-quiet AGNs 
with 0.3 ^ z < 0.96 (upper panel) and with 1.5 ^ z ^ 3.3 (lower 
panel). The values of the X-ray luminosities and spectral indices 
were obtained by fitting the spectra in the observed-frame energy 
range of 0.5-8 keV (see Table \2l . The dashed lines indicate linear 
fits to the data using the least-squares method. The open symbols 
correspond to sources having log A'^h ^ 22, and the filled symbols 
are sources with log Nn > 22. Circles correspond to type 1 AGNs 
and squares to non-type 1 AGNs. 



Qa. 1.70 



(a) 0.5-8 keV Ob served- Frame 


(b) 2-10 keV Resi-Frame 


99.9%... 

68*.- 

0.3 S z S 0.96 

1.5 ^ z ^ 3.3 


0.3 S z S 0.96 

1.5 ^ z ^ 3.3 



0.10 0.15 0.20 0.25 0.30 

a 



0.0 0.1 0.2 0.3 

a 



Fig. 11. — 68% and 99.9% confidence contours of the slope a 
and offset /3* of the P — Lx correlation for AGNs in the 0.3 < 
z < 0.96 (solid contours) and 1.5 < z < 3.3 (dotted contours) 
redshift ranges. The parameters a and /3* were derived from fits 

" log TTn^'"' 



of the linear model F 



(3*. The fits were 



(10** crgs/s) 

performed in the 0.5—8 keV observed-frame band (a) and in the 
2-10 keV rest-frame band (b). The confidence contours indicate 
that the parameters of the linear fit to the F — Lx correlation differ 
at the > 99.9% confidence level between the 0.3 < z < 0.96 and 
1.5 < z < 3.3 redshift ranges. 

0.5-8 keV (see Table [5]). We searched for a correlation 
between F and Lx by computing the Spearman's and 
Kendall's correlations (see Table[3]). We find a strong cor- 
relation between F and L2-10, at the >99.9% confidence, 
for sources having 0.3 ^ z < 0.96 and 1.5 < z < 3.3. We 
tested for a linear dependence between F and log Lx by 
calculating the Pearson's correlation and find a high sig- 
nificance (>99.8%) for sources within 0.3 < z < 0.96 and 
1-5 ^ z < 3.3 (see Table[3]). In Table[Sl we also present 
results of linear least-squares fits to the F — Lx relation 
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TABLE 4 

Correlation table of T versus Lx for AGNs with spectroscopic redshifts. 



r vs. L2-10 

Fitted Energy Frame Redshift bin A^^ rc^ % sign'^ fraction of type 1 fraction with logA^H ^22 



0.5-; 


i keV observed-frame 


0.3 <z <0.96 


46 


0.47 


99.9 


0.37 


0.70 


2- 


-10 keV rest-frame 


0.3 <z <0.96 


40 


0.64 


>99.9 


0.40 


0.68 


0.5-i 


i keV observed-frame 


0.96 < 2 < 1.5 


31 


0.40 


97.3 


0.42 


0.71 


2- 


-10 keV rest-frame 


0.96 < z < 1.5 


26 


0.42 


96.6 


0.50 


0.73 


0.5-! 


i keV observed-frame 


1.5 <z <3.3 


26 


0.38 


94.3 


0.62 


0.62 


2- 


-10 keV rest-frame 


1.5 < z < 3.3 


24 


0.49 


98.5 


0.67 


0.67 



^ Number of RQ AGNs in each redshift bin. Spearman correlation coefficient. ^ Percentile significance of the correlation. 



2.5 - 0.3<z<0.B6 



2.0 ■ 



1 .5 - 




type 1 AGN 
non-type 1 AGN 



42.5 43.0 43.5 44.0 44.5 

log L2_,o 



2.5 
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^3.5 
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Fig. 12. — r versus 2-10 keV luminosity of radio-quiet AGNs 
with 0.3 ^ 2 ^ 0.96 (upper panel) and with 1.5 ^ z ^ 3.3 (lower 
panel). The values of the X-ray luminosities and spectral indices 
were obtained by fitting the spectra in the rest-frame energy range 
of 2-10 keV (see Table [2J. The dashed lines indicate linear fits 
to the data using the least-squares method. The open symbols 
correspond to sources having log A'^e ^ 22, and the filled symbols 
are sources with log Nn > 22. Circles correspond to type 1 AGNs 
and squares to non-type 1 AGNs. 



with a model of the form T = a log Lx + /3. For this 
test, we assumed that F is the dependent variable with 
errors given at the 68% confidence level. In Table El we 
show the best-fit linear fit parameters a and /3. We find 
that the best-fit parameters a and /3 show a significant 
change between the redshift bin 0.3 ^ z < 0.96 and the 
redshift bins of 0.96 < ^ < 1.5 and 1.5 < z < 3.3. In 
particular, for spectral fits performed in the 0.5 - 8 keV 
observed-frame band we find the following: The slope 
and offset of the linear fit to the F — Lx correlation in 
the 0.3 < z < 0.96 redshift range are, a= 0.14 ± 0.02 
and f3 — —4.5 ± 0.8, respectively. The slope and offset 
of the F — Lx correlation in the 1.5 < z < 0.33 redshift 
range are, a — 0.23 ± 0.03, (3 ~ — 8.7± 1.2, respectively. 
Similar result are found for spectral fits performed in the 
2-10 keV rest-frame. This change in the linear parame- 
ters can also be seen in Figure HI 

In Figure [TT] we show the 68% and 99.9% confidence 
contours of a and /3* for AGNs in the 0.3 < 2; < 0.96 



and 1.5 < z < 0.33 redshift ranges and for fits per- 
formed in the 2-10 keV rest-frame band (Figure [TTk ) 
and in the 0.5-8 keV observed- frame band (Figure [TTb). 
The parameter (3* is obtained from fits of the model 

-^2-10 



+ /?*. L 



2-10 



was re-normalized 



T — a log TTTOT J-T 

(10^^ crgs/s) 

for the purpose of illustrating better the full range of the 
contours. The 68% and 99.9% confidence contours levels 
correspond to Ax^(a, (3*) values of 2.3 and 13.81, respec- 
tively. The confidence contours indicate that the param- 
eters of the linear fit to the F — Lx correlation change at 
the > 99.9% confidence level between the 0.3 < z < 0.96 
and 1.5 < z < 0.33 redshift ranges. 

To test the sensitivity and stability of these confidence 
contours to possible outliers in the data we repeated 
the confidence contour analysis by excluding data points 
with significant deviations from the linear fit. In partic- 
ular, we re-fit the F — Lx correlation and re-calculated 
the confidence contours after excluding data points that 
deviated by more than 2a, 2.5a and 3a from the linear 
fit. In all cases we find that the parameters of the lin- 
ear fit to the F — Lx correlation change between redshift 
bins 1 and 3 at the > 99.9% and > 98% confidence levels 
for fits performed in the 0.5-8 keV observed-frame and 
2-10 keV rest-frame, respectively. 

As discussed in §4 to test the influence of possible ef- 
fects such as Compton reflection, soft excesses, and in- 
trinsic absorption on the F - Lx correlation, we also fit- 
ted the spectra in the 2-10 keV rest-frame, where these 
effects are expected to be smaller. The results of these 
spectral fits are presented in Table [21 In Figure [TH we 
present F versus L2-10 for sources in the redshift range 
of 0.3 ^ z < 0.96 (upper panel), and in the redshift range 
of 1.5 ^ z < 3.3 (lower panel) for spectral fits performed 
in the 2-10 keV rest-frame band. The results of our cor- 
relation analysis applied to the variables F and Lx are 
shown in Table [31 We find the Spearman, Kendall and 
Pearson correlation coefficients of F vs. L2-10 to be sig- 
nificant at the >99.9% and >99.7% confidence levels, for 
sources within 0.3 < z < 0.96 and 1-5 < z < 3.3 respec- 
tively. These results suggest that Compton reflection, 
soft excesses, and intrinsic absorption are most likely not 
driving the observed correlation between F and Lx in the 
two redshift bins analyzed in this section. In §5.2.3 and 
§5.3.6, we provide detailed analyses to show that intrinsic 
absorption and Compton reflection have negligible con- 
tributions to the F — Lx correlation. 

5.2.3. Dependence of the F — Lx correlation on Nh 

The estimated values of the photon indices used in 
our correlation analysis depend partially on the assumed 
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TABLE 5 

Results of linear fits to the F vs. log L2-10 relation. 



Redshift bin 


Fitted energy range 


Sample'' 


a 


/3 


0.3 <z <0.96 


observed-frame 0.5-8 keV 


CDFs 


0.14 ±0.02 


-4.5 ±0.8 


0.3 <2 <0.96 


rest-frame 2-10 keV 


CDFs 


0.13 ±0.04 


-4.1 ± 1.7 


0.96 <z < 1.5 


observed-frame 0.5-8 keV 


CDFs 


0.23 ±0.03 


-8.3 ± 1.5 


0.96 %z%1.5 


rest-frame 2-10 keV 


CDFs 


0.27 ±0.05 


-9.9 ±2.2 


1.5 <z <3.3 


observed-frame 0.5-8 keV 


CDFs 


0.23 ±0.03 


-8.7 ± 1.2 


1.5 %z%3.3 


rest-frame 2-10 keV 


CDFs 


0.24 ±0.06 


-8.9 ±2.4 


1.5 < 2 < 3.3 


observed-frame 0.5-8 keV 


combined 


0.27 ±0.03 


-10.3 ± 1.2 



^ Based on fits of a linear model (F — a log Z/2~io + /9) to the F versus log L2-10 relation, 
using the weighted least-squares method. The errors in F at the 68% level are used in the 
linear fit. L2 — 10 is the 2—10 keV rest-frame luminosity. ^ The CDFs sample consists of all 
the sources presented in Table 2. The combined sample consists of sources obtained from 
the independent surveys of lVignali et al. (1999), Geor ge" et al.| ( |2000D , iReeves & Turner I 
j2000D . fPage et al.| j 20031) and [DaTeOQ ll^OOl) combined with the sources of our CDFs 



sample. 



2.0 : 
1 .8 - 

A 

V 1,6 ; 

1 .4 : 

1 .2 _ 
23.5 r 

23.0 r ^19' z 

V 22.5 r z 
- 22.0 r .■t''' z 

21.5 7.-1-''' ^ ^ ^ 

21.5 22.0 22.5 23,0 23,5 
log N„ 

Fig. 13. — Results from fits performed on 1000 simulated spectra 
with S = 550, r = 1.6, 2=1.4 and 5 different values of log A'jj 
(21.5, 22, 22.5, 23 and 23.5). The estimated parameters with fits 
performed in the 0.5—8 keV observed-frame band are shown with 
filled squares (solid error bars), and the estimated parameters with 
fits performed in the 2—10 keV rest- frame band are shown with open 
circles (dashed error bars). A small shift in the value of log A'^h 
(horizontal axis) has been introduced for visual purposes. In the 
upper panel we show (F) vs. log N^. In the lower panel we show 
log(AfH) vs log Afjj. AH error bars represent ± 1 cr deviations. 



TABLE 6 
Results of simulations to test the 
probability of detecting intrinsic 

absorption through SPECTRAL FITS. ^ 



log A^H 


Percentage of cases'" 


21.5 


41.8 % 


22.0 


74.7 % 


22.5 


99.1 % 


23.0 


100.0 % 


23.5 


99.9 % 



^ For each of the 5 different values of log A^h 
(log Nh = 21.5, 22, 22.5, 23 and 23.5), we randomly 
generated 1000 fake spectra assuming sources close 
to the aim-point of the Chandra ACIS-I CCD, with 
550 counts in the 0.5-8 keV band, T — 1.6, and 
z — 1.4. ^ Percentage of eases out of the 1000 sim- 
ulated cases where the i^-test indicates a significant 
presence of absorption assuming the input values of 
column densities listed in the first column. 



spectral models used to fit the AGN spectra. In par- 
ticular, the default model used in our spectral analysis 
assumes a simple power law that can be modified by in- 
trinsic absorption. There is some evidence suggesting 
that the intrinsic column density (N} i) could be evolv- 
ing both with X-r ay luminosity (e.g . lUeda et al.l 20031 ; 



Akvlas et al."2006') and redshift (e.g.. lAkvlas et al.ll2006l : 

Treister & Urry 2006). At the same time, large values 
of A^H could be producing some dispersion in the esti- 
mated values of F. In order to analyze the effect of A^h in 
the spectral fitting, we have performed simulations using 
the software command "fakeit" in XSPEC. We randomly 
generated 1000 fake spectra for each of five values of A'h 
(log Nil = 21.5, 22, 22.5, 23 and 23.5). Each simulated 
spectrum was created assuming an absorbed power-law 
(APL) model with 550 counts in the 0.5-8 keV band, 
r — 1.6, and z — 1.4. The simulated sources were con- 
sidered close to the aim-point of the Chandra ACIS-I 
CCD. The assumed values of the total counts, T and z are 
close to the mean values found in §4. We performed fits 
to the randomly generated spectra using the same APL 
model, and plot in Figure [T3] (upper panel) , the mean 
spectral slope (with standard deviation) of the 1000 fits 
as a function of log A'h; these fits were performed both 
in the 0.5-8 keV observed-frame (squares) and the 2- 
10 keV rest-frame (circles). Based on these results we do 
not find any significant bias in the estimation of F with 
A'h- We do, however, find that the standard deviation 
shows a clear tendency to grow with Ajj independently of 
the energy band fit, as seen in Figure fT3l (upper panel). 
In Figure [13] (lower panel), we see that in general the 
estimated value of A^h is accurate for log Ah > 22; how- 
ever, for log A^H < 22 the column density is slightly 
overestimated and has a larger dispersion. 

Using the same simulations, we estimated the effective- 
ness of using the i^-test at the 95% level of significance 
to determine the improvement in the fit quality by using 
an absorbed power-law (APL) model as an alternative to 
the default power-law (PL) model. Table[S]shows that in 
a simulation of 1000 fake spectra with S — 550, F = 1.6, 
z — 1.4 and log A^h = 22, the i^-test indicates absorp- 
tion in ~75% of the spectra. For simulated spectra with 
log Ah = 22.5, the i^-test indicates absorption in ^ 99% 
of the cases and for log A'h > 22.5 the F-test indicates 
absorption in more than 99.9% of the cases. Based on 
these simulations, we conclude that the F-test can accu- 
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TABLE 7 

Correlation table of F vs. Lx for sub-samples of different absorption. 



Redshift bin 


Absorption" 


Fitted energy range 


A*" 


L2~ 
re" 


10 vs. F 
% sign'^ 


n S < 7 < n Qfi 


log 


JVj] < 22.5 


0.5- 


3 keV observed-frame 


40 


0.38 


98.4 


3 < 7 < n Qfi 


log 


JVj] < 22.0 


0.5- 


3 keV observed-frame 


37 


0.38 


98.1 


3 < 7 < n Qfi 


log 


JVj] > 22.0 


0.5- 


3 keV observed-frame 


16 


0.66 


99.4 


s < ? < n Qfi 






0.5- 


3 keV observed-frame 


17 


0.34 


81.6 


n s < ? < Qfi 


log 


JVh < 22.5 


2 


— 10 keV rest-frame 


32 


0.56 


>99.9 


3 < z < Qfi 


lop" 
log 


Afjj 22.0 


2 


—10 keV rest-frame 


27 


0.56 


99.8 


n 3 < 7 < Qfi 


log 


jVt, > 22 


2 


— 1 rt I'fiQt— Tm m f> 
1-KJ jYt; V 1 1;.!) L 11 diiic; 


17 


0.67 


99.7 


3 < 2 < n 96 




2 


— 10 keV rest-frame 


16 


0.58 


98.2 


Qfi < 2 < 1 5 


lop" 

lug 


Nil ^ 22.5 


0.5- 


3 keV observed-frame 


38 


0.27 


89.9 


n Qfi < r < 1 ^ 


log 


Nil < 22.0 


0.5- 


8 keV observed-frame 


31 


0.31 


91.6 


n Qfi < r < 1 


log 


TVh > 22.0 


0.5- 


8 keV observed-frame 


23 


0.31 


85.4 


Qfi < 2 < 1 5 




type 1 


0.5- 


8 keV observed-frame 


17 


0.20 


49.5 


0.96 <2 < 1.5 


log 


A^H < 22.5 


2 


-10 keV rest-frame 


33 


0.26 


85.3 


0.96 <2 < 1.5 


log 


Ah < 22.0 


2 


-10 keV rest-frame 


27 


0.33 


91.1 


0.96 < 2 < 1.5 


log 


Ah > 22.0 


2 


-10 keV rest-frame 


19 


-0.18 


54.8 


0.96 < 2 < 1.5 




type 1 


2 


-10 keV rest-frame 


16 


0.15 


37.1 


1.5 <z £3.3 


log 


Ah < 22.5 


0.5- 


8 keV observed-frame 


40 


0.45 


99.7 


1.5 <2 <3.3 


log 


Ah < 22.0 


0.5- 


8 keV observed-frame 


32 


0.35 


95.2 


1.5 <2 <3.3 


log 


Ah > 22.0 


0.5- 


8 keV observed-frame 


25 


0.38 


94.2 


1.5 <2 <3.3 




type 1 


0.5- 


8 keV observed-frame 


16 


0.12 


34.0 


1.5 <2 <3.3 


log 


Ah < 22.5 


2 


-10 keV rest-frame 


34 


0.45 


99.2 


1.5 ^2 ^3.3 


log 


Ah % 22.0 


2 


-10 keV rest-frame 


23 


0.42 


95.6 


1.5 <2 <3.3 


log 


Ah > 22.0 


2 


-10 keV rest-frame 


23 


0.42 


96.2 


1.5 <2 <3.3 




type 1 


2 


-10 keV rest-frame 


16 


0.29 


72.1 



^ Sub-samples contain RQ AGNs that arc cither type 1 AGN or have Nh less than a 
specified value. ^ Number of RQ AGNs in eaeh sub-sample. The Spearman correlation 
coefficient.^ The significance of the Spearman correlation coefficient. 



rately identify absorption when log Nn > 22. 

Based on our finding that highly-absorbed sources 
show a greater dispersion of the estimated value of T, 
we tested the sensitivity of the F — Lx correlation for 
sources having 0.3 ^ z < 0.95 and 1-5 < z < 3.3 to in- 
trinsic absorption, by removing sources with significant 
absorption (log A'h ^ 22.5). We also tested this cor- 
relation for sources having log A^h < 22. Finally as a 
complementary test we analyzed the F — Lx correlation 
for type 1 AGNs. The results of these three tests are pre- 
sented in Table [71 For sources having log A'h ^22.5 and 
log A'h ^ 22 we find in the first and third redshift bins 
that the Spearman correlation coefficients of F versus 
-^2-10 are significant at the >95% confidence levels. This 
result holds for sources with fits performed in both the 
0.5-8 keV observed-frame and in the 2-10 keV rest-frame 
(see Table III). Notice that sources having log iVn ^22 
are plotted as empty squares in Figures \W\ and [121 

For type 1 AGNs, we find that the F— L2-10 correlation 
is significant at the 82% and 12% levels in the first and 
third redshift bins, respectively, for fits performed in the 
0.5-8 keV observed-frame band; the significances are at 
the 98% and 72% levels, respectively, for fits performed 
in the 2-10 keV rest-frame band. We briefly investigate 
possible reasons that may explain the apparent low de- 
tection significance of the F — L2-10 relation for the type 
1 AGNs of our sample. First we note that the lumi- 
nosity ranges of the type 1 AGNs of our sample in the 
first and third redshift bins are 42.6 < log £2-10 S 44.7 
and 43.8 < logL2-io ^ 44.8. Our sample of type 1 
AGNs therefore includes relatively luminous sources in 
each redshift bin. For sources in the third redshift bin, 
as we will later see in §5.2.7, the values of F appear to 
saturate above log L2-10 ~ 45. Therefore, the type 1 



2.5 - 




1.6 1.8 2.0 2.2 2.1 2.6 2.8 3.0 
2 

Fig. 14. — Spectral slope (F) versus redshift (2) of radio-quiet 
AGNs with 1.5 < 2 < 3.3 with fits performed in the 2-10 keV 
rest-frame. The size of the symbols increases with L2-10- 

AGNs detected in the third redshift bin of our sample 
are expected to lie on the fiat part of the F — Lx re- 
lation. We conclude that the apparent low significance 
of the F — ^2-10 relation for the type 1 AGNs of our 
sample found in the third redshift bin is mainly the re- 
sult of their relatively large luminosity and the limited 
number of type 1 AGNs in our sample. Based on the 
tests presented in this section we confirm that the strong 
F — Lx correlations that we find in RQ AGNs in the red- 
shift ranges of 0.3 < 2: < 0.95 and 1-5 < z < 3.3 are not 
driven by intrinsic absorption. 

5.2.4. Other Correlations 

We searched in the redshift ranges of 0.3 < z < 0.95 
and 1.5 < z < 3.3 for possible correlations between the 
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TABLE 8 



Correlation table of L2- 


-10 vs. z, r vs. 


2, A'"h vs. z, Na vs. -L2-10 


AND r vs 




correlated parameters 




Redshift bin 


Fitted energy Range 




% sign*^ 


L2-10 vs z 


53 


{j.o ^ z ^ U.yD 


U.o— keV observed- irame 




yo.o 


r VS 2 


53 


{j.o ^ z ^ u.yo 


0.5-^8 keV observed- frame 


U. LZ 


f;i /I 

01 .4 


Afjj vs z 


27 


U.o ^ 5; ^ U.yo 


0.5-^8 keV observed- frame 


— U.Uo 


Qn 
ou.y 


iVn vs L2-10 


27 


n ^ ^ n Qfi 
u.o ^ ^ ^ u.yo 


0.5—8 keV observed-frame 


u. 10 


ou. / 


r vs A^H 


27 


U.o ^ -0 ^ u.yo 


0.5—8 keV observed- frame 


n 1 /I 


E^Q « 

00.0 


-^^2-10 VS Z 


44 


n <^ < n Qf; 
U.vS ^ -2 ^ u.yo 


2—10 keV rest-frame 


n 00 




r vs 2 


44 


\j.o ^ z ^ U.yo 


2—10 keV rest-frame 


U.Zo 


fi7 /I 
01.4 


A^H vs 2 


25 


n ^ <" n Of; 
U.o ^ -2 ^ U.yo 


2—10 keV rest-frame 


— U.Uo 


on f; 

ZU.O 


A^H VS L2-10 


25 


n ^ n Of; 
U.o ^ -2 ^ U.yo 


2—10 keV rest-frame 


U.Uo 


Qn 
oU.z 


r vs A^H 


25 


u.o ^ z ^ u.yo 


2 — 10 keV rest-frame 


u.oz 


00.0 


^2-10 vs 2 


54 


n Of; ^ -» <r 1 
U.yo ^ -3 ^ i.o 


0.5—8 keV observed- frame 


n 1 f; 

U. ID 


7E^ f; 
( 0.0 


r vs 2 


54 


n Of; ^ -V <r 1 
U.yo ^ 2 1.0 


0.5—8 keV observed- frame 




yi.o 


Afjj vs 2 


28 


n Of; ^ -V <r 1 E^ 
U.yo ^ 2 i.o 


0.5—8 keV observed- frame 


— U.oo 


01 c 
yi.o 


A^H vs L2_10 


28 


n Of; ^ -V <r 1 E^ 
U.yo ^ z 1.0 


0.5—8 keV observed- frame 


— U. lo 


00.4 


r vs A^H 


28 


n Of; ^ -V <r 1 E^ 
U.yo ^ -S ^ 1.0 


0.5—8 keV observed- frame 


— U.oo 


00 7 
yz. i 


L2-IO vs 2 


46 


n Of; ^ ? <r 1 E^ 
u.yo ^ -3 ^ -L-O 


2—10 keV rest-frame 


u. 1 / 


7A (\ 
/ 4.0 


r vs 2 


46 


n Of; ^ -V <r 1 E^ 
U.yo ^ 2 i.o 


2—10 keV rest-frame 


— U.Uo 


/II 1 
41.1 


Afjj vs 2 


24 


n Of; ^ -» <r 1 e:; 
U.yo ^ -3 ^ 1.0 


2—10 keV rest-frame 


— U.oD 


01 Q 

yi.o 


A^H vs L2-10 


24 


n Of; ^ -» <r 1 e:; 
U.yo ^ -3 ^ 1.0 


2—10 keV rest-frame 


— U. iU 


Qf; n 

oD.U 


i vs iVH 


24 


n Of; ^ -» <r 1 
U.yo ^ -3 ^ I.o 


2—10 keV rest-frame 


— U.UvS 


11.0 


L2-10 vs z 


57 


1.5 < 2 < 3.3 


0.5-8 keV observed-frame 


0.58 


>99.9 


T vs z 


57 


1 'i ^ 7 <! 3 3 


W.O IVC. V ULJo^i V d_i il. (Xlil^ 


0.14 


68.4 


TVh vs z 


34 


1.5 < ^ < 3.3 


0.5-8 keV observed-frame 


0.13 


53.5 


Nil vs L2-10 


34 


1.5% z% 3.3 


0.5-8 keV observed-frame 


-0.23 


81.8 


r vs A^H 


34 


1.5 < ^ < 3.3 


0.5-8 keV observed-frame 


-0.11 


48.3 


^'2-10 vs Z 


48 


1.5 < ^ < 3.3 


2—10 keV rest-frame 


0.65 


>99.9 


r vs 2 


48 


1.5 < ^ < 3.3 


2—10 koV rest-frame 


0.35 


98.6 


A^H vs z 


30 


1.5% z% 3.3 


2-10 kcV rest-frame 


0.34 


93.4 


Nb. vs L2-10 


30 


1.5 < ^ < 3.3 


2—10 keV rest-frame 


-0.02 


8.8 


r VS Nil 


30 


1.5% z% 3.3 


2-10 keV rest-frame 


-0.1 


40.4 



^ Number of RQ AGNs in each sub-sample. The Spearman correlation coefficient.'^ The significance 
of the Spearman correlation coefficient. 



photon index T and other physical parameters of our 
sample of AGNs by computing Spearman rank corre- 
lations. The results from this correlation analysis are 
shown in Table[8l In the cases where the intrinsic column 
densities were TVh ^ 10^° cm~^, we could only obtain up- 
per limits to A'h, and therefore computed the correlation 
coefficients using survival analysis (Isobe ct al. 1986). 

The selection criteria used in this work impose a lumi- 
nosity limit which is redshift dependent (see Figure [5]) . 
Furthermore, the co-moving density of luminous AGNs is 
known to increase with z. These two effects will produce 
a correlation between Lx and z as can be seen in Table [51 
We do not find any significant correlation between F and 
A'h in any bin. The fact that F is not correlated with 
Nh provides further support that the F — Lx correlation 
is not driven by Njj- 

We also find a weak correlation between A^h and z in 
the redshift range of 1-5 < z < 3.3 for fits in the 2- 
10 keV rest-frame band. This result may imply that the 
intrinsic column density evolves, increasing with redshift. 
Such a result h as been reported in several studies (e.g., 
iLa Franca et a l. 2005; Trcister & Urry 2006); however, 
the evolution of A^h with z is still a debatable topic since 
other authors have not found definitive evid ence for the 
evolution in the "obscura.tion fraction" (e .g..lUeda et al.l 
[200I lAkvlas et al.ll200llDwellv k Page \\2Qm . 

A rather surprising result was the detection of a corre- 
lation between F and z in the third redshift bin with an 
apparent significance of 99.6%. A careful analysis indi- 
cates that this apparent T — z correlation is most likely 
the result of selection effects. This tendency seems to be 



confirmed in Figure [TH This F versus z plot indicates 
that higher luminosity sources tend to group in the upper 
right area and lower luminosity sources in the lower left 
area. To test for selection effects, we performed a correla- 
tion analysis including sources with luminosities greater 
than the minimum luminosity of a detectable source at 
z 3. This limit corresponds to log ^2-10 ^ 44 (see Fig- 
ure[5]). We find that the Spearman's correlation probabil- 
ity of the T — z correlation in the third redshift bin for fits 
performed in the 2-10 keV rest-frame, decreases to a non- 
significant level of ~ 66% when we only include sources 
with log L2-10 > 44 (30 RQ AGNs). We note, however, 
that within the same luminosity range the F — L2-10 cor- 
relation is significant at the >99% {rc ^ 0.6) level. Our 
analysis indicates that the apparent correlation between 
F and z in the third redshift bin is most likely the result 
of selection effects. This conclusion is also confirmed in 
§5.2.7. 

5.2.5. Radio-Quiet AGNs with 0.96 < z < 1.5 

In §5.2.1, we showed that the F — Lx correlation was 
not significant for sources having 0.96 ^ z < 1.5, espe- 
cially for fits performed in the 2-10 keV rest-frame where 
we found that the correlation was only significant at the 
~85% confidence level. In this section we investigate the 
cause of the lower significance of the F — Lx relation for 
sources having 0.96 ^ z < 1.5. In Figure [T5l we show 
F versus ^2-10 for sources having 0.96 ^ z < 1.5 with 
fits performed in the 0.5-8 keV observed-frame (upper 
panel) and 2-10 keV rest-frame (lower panel). The F ver- 
sus £2-10 data points show a larger scatter than what is 
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Fig. 15. — F versus 2—10 keV luminosity of radio-quiet 
AGNs having 0.96 ^ z ^ 1.5 with fits performed in the 0.5—8 keV 
observed-frame band (upper panel) and fits performed in the 
2-10 keV rest-frame band (lower panel). The dashed lines indi- 
cate linear fits to the data using the least-squares method. The 
open symbols represent sources having log Nn < 22, and the filled 
symbols represent sources with log Afjj > 22. Circles represent type 
1 AGNs and squares non-type 1 AGNs. 

seen in the other redshift bins consistent with the lower 
significance found for the F — Lx correlation. In Ta- 
ble [71 we present the results of our correlation analysis 
of the r — Lx data for sub-samples of different intrin- 
sic absorption. We find that sources with log A'h 22 
and 0.96 ^ ^: ^ 1-5 show no significant correlation be- 
tween r and Lx, whereas sources with log Nn < 22 and 
0.96 < 2 < 1.5 have a F — Lx correlation that is signifi- 
cant at the >90% confidence level. We conclude that the 
absorbed sources with 0.96 ^ ^ < 1-5 are possibly dilut- 
ing the correlation significance found in this redshift bin. 
We caution, however, that the low number of sources per 
sub-sample used in this analysis combined with the un- 
certainties in the photon indices may result in inaccurate 
estimates of the strengths of the correlations. 

Our correlation analysis between several other spec- 
tral parameters for sources with 0.96 ^ ^ < 1-5 is in- 
cluded in Table [51 We find an anticorrelation between 
A'h and z at the ^ 90% confidence level for sources with 
0.96 ^ z < 1.5 and for fits performed in the 0.5-8 keV 
observed-frame and the 2-10 keV rest-frame bands. Cor- 
relations between T versus z and T versus iVu for sources 
in 0.96 ^ z <1.5 are found to be moderately significant 
for fits performed in the 0.5-8 keV observed-frame band 
and not significant for fits performed in the 2-10 keV 
rest-frame band. 

5.2.6. Dependence of the F — Lx correlation on 
Compton-re flection 

In this section, we address the possibility that the 
r — Lx correlation found in this work is produced by a 
change with luminosity of the Compton-reflection com- 
ponent. We note that the Compton-reflection compo- 
nent is difficult to model accurately in low-to-medium 
S/N X-ray spectra and therefore inaccurate modeling of 
this component may result in apparent fiattening of the 
X-ray spectra. 

Several studies indicate that the equivalent width 



(EW) of the iron Ka emission line in the X-ray spec- 
tra of AGNs is anti-correlated with the 2-10 keV lumi- 
nosity (e.g.. jlw asawa fc Taniguchi ' '1993l: _|Nandra et al.l 
Il997t iPage etal . 2004a; Bianchi et al. 2003). This anti- 
correlation is commonly referred to as the 'X-ray Baldwin 
effect' and is also known as the Twasawa & Taniguchi 
effect' . There are several proposed physical explana- 
tions in the literature for the X-ray Baldwin effect in- 
cluding (1) a change in the coveri ng factor of a Compton - 
thick torus with luminosity fe.g.. lKomgl fc Kartielll994[) . 
(2) a luminosity-depen dent ionization stat e of the iron- 
emitting material (e.g.. lNandra et afll997t l. and (3) vari- 
ability of the continuum AGN emission as suming con- 
stant iron-fine fluxes (e.g.. lJiang et al.|[2006[ ). It has also 
been proposed that the X-ray Baldwin effect is driven 
mostly by changes in the Eddingto n luminosity ratio 
rather than by X-ray luminosity fe.g. JJiang et al.| [2006). 

Several models of AGN accretion disks assume the iron 
line and Compton-reflection components originate from 
X-ray emission reprocessed in the accretion disk and indi- 
cate that the strength of the Compton-reflection compo- 
nent increases monotonically with the EW of the iro n line 
(;e.g.. lGeorge fc Fabian l[l99ll : iGhisellini et al.lll994f) . and 
consequently decreases with Lx as well. Therefore, un- 
der this premise, these models could possibly explain the 
F — Lx relation found in this work, since a decrease of the 
Compton-reflection component with Lx could result in 
an increase in F with Lx if the Compton-reflection com- 
ponent is not modeled accurately in our spectral analysis. 

Bianchi et al. (2007) recently found a strong anti- 
correlation between the neutral narrow component of 
the iron Ka emission line and the 2-10 keV luminosity 
of AGNs. These authors suggest that the neutral nar- 
row iron-line component originates from the molecular 
torus and the broad iron-line component originates from 
reprocessing in the accretion disk. The dependencies, 
however, of the Compton-reflection compo nent with the 
2-10 keV luminosity are still unclear fe.g.. lNandra et ahl 
Il995t iPage et al.ll2004bl ). It is also unclear whether the 
Compton-reflection component observed in the X-ray 
spectra of AGNs originates mostly from the torus or the 
accretion disk. If the neutral narrow iron Ka emission 
line originates from the torus and the Compton-reflection 
component from the accretion disk then one cannot sim- 
ply assume that the Compton-reflection component will 
follow the X-ray Baldwin effect. Variability studies of 
indivi dual AGNs such a s the Seyfert 1 gala xies NGC 
5548 (iChiang et al.l [20M) and MCG-6-30-15 (|Lee et all 
I2OOOI) indicate that the Compton-reflection component 
increases with X-ray luminosity and the iron line EW 
and the relative normalization of the Compton-reflection 
hump are anti-correlated. We note that the Seyfert 1 
galaxies in these variability studies were observed to vary 
over a factor of up to w3 in luminosity whereas our study 
includes objects spanning a factor of ~200 in luminosity. 
It is therefore difficult to extrapolate the results of these 
variability studies to our work. 

Since observationally it is still unclear how the 
Compton-reflection component depends on X-ray lu- 
minosity we have investigated the degree to which 
Compton-reflection can be driving the F — Lx relation 
by performing simulations and additional tests upon 
our data. We first simulated X-ray spectra containing 
Compton-reflection components with integer reflection 
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Fig. 16. — Results from fits performed on 1000 simulated spec- 
tra with 5=550, r=1.9, and 5 different values of the Compton- 
reflection scaling factor R (il=0,l,2,3,4). The mean of the best-fit 
values of F with fits performed in the 0.5—8 keV observed-frame 
band are shown with filled squares (solid error bars) , and the mean 
of the best-fit values of F with fits performed in the 2—10 keV rest- 
frame band are shown with open circles (dashed error bars). A 
small shift in the values of R (horizontal axis) has been introduced 
for visual purposes. In the upper panel we show (F) vs. R with 
2=0.7. In the lower panel we show (F) vs. R with z=2.2. On 
the X-axis we also show the ratio of photons (/ij) in the full band 
(0.5-8 keV) that originate from Compton reflection to photons that 
originate from the direct power-law component. All error bars rep- 
resent ± 1 (T deviations. 

scaling factors ranging between < i? < 4. For each 
value of the reflection scaling factor (i?=0,l,2,3,4); we 
simulated 1000 spectra using the FAKEIT command in 
XSPEC. The Compton-reflection components were sim- 
ulated using the PEXRAV model, assuming sources close 
to the aimpoint of the ACIS-I CCD, r=1.9, log TVr = 22, 
a total number of events per spectrum of S'=550, an e- 
folding cutoff energy of i?cut=400 keV, and an inclination 
angle (i) of the reflector equal to 30°. The values of Ecut 
and i were chosen to be close to those generally used 
to m odel Seyfert galaxies (e.g.. iMagdziarz fc Zdziarskil 
119951 ). Our results are insensitive to any reasonable 
value of Ecut and changing i will mostly affect the over- 
all strength of the reflection component. We performed 
these simulations assuming redshifts of z — 0.7 and z = 
2.2, which correspond to the mean redshifts of the sources 
in our sample with 0.3 ^ z < 0.96 and 1-5 < z < 3.3. We 
proceeded in fitting the simulated spectra with absorbed 
power-law models to estimate the decrease in the fitted 
values of F vs. the strength of the Compton-reflection 
component. In Figure [TBI we show the best-fit values of 
F as a function of the reflection scaling factor for sources 
with redshifts of z = 0.7 (upper panel) and z — 2.2 (lower 
panel). In Figure [TBI we also show the ratio of photons in 
the full band (0.5-8 keV) that originate from Compton 
reflection to photons from the direct power-law compo- 
nent (/r). 

As expected fits performed in the 2-10 keV rest-frame 
band of the high-z sources are less affected by the 
Compton-reflection component and show a smaller 
change of F than fits performed in the 0.5-8 keV 
observed-frame band. Specifically, we find apparent 
changes of F of about 0.7 and 0.3 for fits performed in 
the 0.5-8 keV observed-frame and 2-10 keV rest-frame 
bands, respectively, for sources with 1.5 ^ z < 3.3. 



Our simulations indicate that if Compton reflection 
is producing the observed change in F of about 0.5 for 
sources with 1-5 < z < 3.3 (see lower panel of Figure [TP)) , 
then the mean values of F derived from fits performed 
in the 0.5-8 keV observed- frame band should differ by 
about 0.2 from the mean values of F derived from fits 
performed in the 2-10 keV rest-frame band. Our obser- 
vations indicate that this is not the case. For sources 
in the redshift bins of 0.3 < z < 0.96, 0.96 < z < 1.5 and 
1-5^2^ 3.3 the differences between the weighted mean 
values of F {{Trest)-{robs)) obtained from fits performed 
in the 0.5-8 keV observed-frame and 2-10 keV rest-frame 
band are -0.03±0.02, 0.02±0.02 and 0.01±0.03 (1-cr er- 
rors), respectively; the similarity between (Trest) and 
(Fobs) is consistent with the results found in §4. Accord- 
ing to our simulations, if the F — ix correlation were pro- 
duced by the Compton-reflection component then these 
differences in the weighted mean values of F would in- 
crease with redshift, reaching values close to 0.2 for 
sources with 1.5 ^ z < 3.3. 

We also expect that if Compton reflection is driving 
the observed F — Lx relation then the strength and slope 
of the correlation for fits performed in the 0.5-8 keV 
observed-frame band should be significantly stronger and 
steeper than the strength and slope of the correlation 
for fits performed in the 2-10 keV rest-frame band, es- 
pecially for the high-redshift sources. This is not the 
case. As shown in Table [31 the strength of the F — Lx 
relation in the 0.5-8 keV observed-frame and 2-10 keV 
rest-frame bands for sources with 1.5 < z < 3.3 is 0.45 
(at the 99.9% confidence level) and 0.43 (at the 99.8% 
confidence level) , respectively. From Table [5] the ratio 
of the slopes of the F — Lx relation for sources with 
1.5 < z < 3.3 is aobs/o^rcst = 0-98 ± 0.25, where a^hs 
and Qfrcst are the slopes derived from fits performed in 
the 0.5-8 keV observed-frame and 2-10 keV rest-frame 
bands, respectively. Based on our simulations the ra- 
tio of the slopes of the F — Lx relation for source with 
1.5 < z < 3.3 would be approximately Ofobs/arcst ^ 1-5 
if Compton-reflection was driving the correlation. 

5.2.7. Radio-Quiet AGNs from other surveys with 
1-5 <« < 3.3 

We also tested the significance of the F — Lx corre- 
lation in the redshift range ofl.5<z<3.3by combin- 
ing results fro m the independ e nt survey s of Vignali et a Ll 
(1999) . iGeorge et all (120001) . [Reeves fc Turner I (|2000f) . 
Page et al.l (|2003[ ). and iD ai et al.[ (|2004[ ) with the re- 
sults obtained from the CDFs. In Figure [T71 we show 
F versus L2_io for sources having 1.5 ^ z < 3.3 com- 
bining the results of these surveys with our fits in the 
0.5-8 keV observed- frame (Table 12). Sources from our 
survey fill in the low-luminosity range of the combined 
data Lz-in ^ lO'*^ - 10"*^ erg s"^ In this range F in- 
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creases with Lx] however, for L2-10 > 10 °erg s 
appears that this relation begins to saturate. 

We find that the Pearson linear correlation coefficient 
reaches a maximum value for sources with L2-W in the 
range of 10'*'^ — 2 x lO**^ erg s~^, and the Spearman rank 
coefficient reaches a maximum value for sources with 
L2-10 in the range of 10^^ — 8 x 10^^ erg s~^. There 
are 76 sources with L2-W in the range of 10'*'^ — 2 x 10^^ 
erg s~^. The Pearson's correlation coefficient for F — Lx 
in this luminosity range is Tp 0.52 and is significant at 
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Fig. 17. — Photon indices of 1.5 < z < 3.3 radio-quiet AGNs 
obtained from several surveys versus their 2-10 keV luminosi- 
ties. Open stars are from Vignali et al. (1999), fi lled stars are 
from'George et al.' f2000), filled circles arc from Reeves fc TurneFI 
lUOOO), open circles arc from Page ct al. (2003), open squares are 
from Dai et al., C2004) and filled squares show data from this work. 
The vertical dashed line corresponds to L2— lO = 2 X 10*^ erg s~^. 
The dot-dashed line shows the best-fit linear model over the lumi- 
nosity range of lO**^ — 2 X 10'*^ org s~^. 

the >99.9% confidence level (null hypothesis probability 
1.8 X 10~^). For reference, the values of the parameters 
of the best-fit model of the form T — a log L2-10 + P can 
be found in Table [S] There are 84 sources with L2-10 in 
the range of 10^'^ — 8 x 10**^ erg s~^. The Spearman cor- 
relation coefficient in this luminosity range is rs ^ 0.60 
and is significant at the >99.9% confidence level (null 
hypothesis probability ^ 1.6 x 10~^). The limits of the 
optimized ranges are marked in Figure 1171 Figure [17] 
suggests two different luminosity regimes of the F — Lx 
relation. In the first regime that covers the luminosity 
range of L2-10 ~ lO''^ — 2 x 10^^ erg we find a linear 
relation between F and log i2-io- In the second regime, 
w here L^-m^ > 2 x 10^^ erg s~^, we confirm the finding 
of iDai et all (|2004D that F decreases with L2-io- Specifi- 
cally, for the 16 sources with L2-10 ^ 2 x 10^^ erg s~^, we 
found F and ix to be anti-correlated with a Spearman 
correlation coefficient of rs ~ —0.50 that is significant 
at the 95.4% confidence level. We also found that for 
the 84 X-ray luminous sources with £2-10 in the range 
of 10'*'^ — 8 X 10^^ erg s~^ F and z are not correlated 
with a Spearman correlation coefficient of ~ 0.09 that is 
significant at the ^60% confidence level. 

5.2.8. Physical Interpretation of the V — Lx Relation and its 
Possible Evolution 

AGN X-ray variability studies have proven very useful 
for improving our understanding of the physical struc- 
tures that produce AGN X-ray spectra. Several vari- 
ability studies of individual AGN have found a positive 
correlation between F and the X-ray lumin osity, (e.g., 
pagdziarz ct al. 1998; Zdziarski "erall |2003[ ) . However, 
the slope of the F — Lx correlation found in the studies 
of individual AGNs appears to be significantly steeper 
than the F v ersus Lx slope f ound in our current study. 
For example, IZdziarski et al.l (|2003( ) report that a change 
in Lx by a factor of 10 results in an increase in F of ~ 
0.6. For a similar change in Lx, we find an increase in F 
of ~ 0.2-0.3 (see Table [5]). 

To explain the F — Lx correlation, we present two 



steady stat e corona models pr o posed by iHaardt et al.l 
(|1997if l and iMerloni fc Fabian I ()2001h . These coronal 
models were originally aimed to explain various X-ray 
variability observations of individual objects and there- 
fore assume a constant black-ho le mass. We also intro- 
duce a third model ( Merloni fc Fabian 1120021 ) focused on 
the relative strength of the corona emission and the ac- 
cretion rate. The main goal for using this model is to 
provide a possible connection for the results obtained 
in our work wi t h the F — Lboi/^Edd correl ations found 
by IWang et al.l (|2004D and IShemmer et"al] p006.) . At 
the end of this section we also comment on how the 
F — Lx correlation depends on the optical depth of the 
hot corona and its evolution with z. 

The first model posits that the inner accretion disk 
is sand wiched by a hot, t enuous and possibly patchy 
corona (|Haardt et al.lll997l ). The corona is coupled to 
a cooler optically-thick layer (accretion disk), which pro- 
vides the seed soft photons that cool the hot layer via 
inverse Compton scattering. The spectrum of the scat- 
tered photons is in general well fitted by a power-law and 
accounts for a large fraction of the observed X-ray emis- 
sion in AGNs. This model predicts that F increases with 
the optical depth of the corona r and it decreases with 
the temperature of the corona. Moreover, if the corona 
is dominated by pairs, the optical depth of the hot 
phase is determined by the compactness I alone ^, where 
t is defined as: 

mgC-^ K r 

Here C is the luminosity of the corona in Eddington 
units, and r is its size in units of Schwarzschild radii. 
For the case of a pair-dominated corona this first model 
predicts that F will increase with Lx- In particular, 
a change in Lx of a factor of 10 is predicted to pro- 
duce a AF 0.2. This predicted variation is consistent 
with that detected in our sample but slightly lower than 
that observed in variability studies of individual AGNs. 
When £ > 0.1 r (^ > 1000) the model indicates that r 
saturates and does not incr ease beyond that luminosity 
(|Haardt fc Maraschi I Il993f) . The first model therefore 
predicts a flattening of the F — Lx relation for bright 
sources as can be seen in Figure [TTl 

The second model originally proposed by 
IMerloni fc FabiaiTI (|2001f ) assumes a static patchy 
corona and is often referred to as the "thundercloud" 
model. The X-ray spectrum in this model is produced 
by thermal Comptonization in spherical regions that 
are raised above the disc at a given height due to 
magnetic flares in active regions of the corona. The 
thundercloud model predicts a F — Lx correlation that 
is consistent with the one we detect for high-redshift 
radio-quiet AGNs in the case where the luminosity of an 
active region increases with its increasing size at a given 
height. Larger active regions tend to be more luminous, 
cooler and produce softer spectra (i.e., F increases). 
If the size of the active region gets too large, then a 
saturation point is reached and the F — Lx relation 
becomes flat for luminous sources in agreement with the 
combined data sample (Figure [T7|). The luminosity of an 

^ pa ir production b ecomes important for ^ -ij 10 (e.g., see 
page 51 of lPetersonHigOTI) 
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active region is assumed to scale with its size r via the 
relation L{r) oc r^, where D \s a. free parameter in the 
model which may be related to the internal structure 
of the region and/or radial dependence of the energy 
generation in the accretion disk. Any change in D will 
directly affect the slope of the F — Lx relation, making 
this model more flexible to explain both our F — Lx 
correlation in high redshift AGNs and those found in 
variability studies of individual AGN (iZdziarski et alj 
[2OOI . 

Under the assumption that the correlation is produced 
by sources of similar black-hole masses, it seems nat- 
ural that changes in X-ray luminosity may result from 
changes in the Eddington ratio e ~ 0.1m, where m is 
the accretion rate in units of AfEdd — -^r*-- Specifi- 
cally, assuming Lx increases with iboi (see equation 21 
of Marconi et al. 2004) and since L\io\ — ^-^Edd 
pect, for sources with similar mass, an increase in e to 
result in an i ncrease in £x- A more d etailed analysis 
is provided in iMerloni &: Fabianl (|2002f ) in the case of 
a coronal-outflow dominated accretion disk model. Un- 
der the assumption that the total power released from 
the accretion disk-corona system is Xboi = e^Edd, a frac- 
tion / will be released in the corona Lx ~ f ■ Lboi- As- 
suming that magnetic turbul ence is the main source o f 
angular momentum transport . IMerloni fc Fabianl ()2002f ) 
concluded that the relationship between / and e can be 
approximately modeled as a power-law / oc with 
S ~ 0.4. This relation is m ostly independent of the 
mass of the black hole Mbh (|Wang et al.ll2004D . As- 
suming that Mbh is kept constant then Lx oc e^~^ , and 
therefore Lx increases with e. Using these relations in 
combination with a steady corona model like the ones al- 
ready describe d in this section, wh i le ass uming that r in- 
creases with e, IMerloni fc Fabianl (j2002f ) concluded that 
the corona gets cooler and F increases with e. This could 
give a plausible explanat ion for the corre l ation between 
F — Lboi/LEdd found bv IShemmer et all ()2006D . Based 
on the above, and assuming similar black-hole masses, 
variations in the Eddington ratio could also explain the 
correlation found in our work since the coronal-outflow 
dominated accretion disk model predicts that F increases 
with Lx- 

An alternative explanation is that the F — Lx corre- 
lation is driven by variations in the black-hole masses 
of the sources and e is constant, which is in agree- 
ment with the predictions of semi-analytic models by 
iKauffmann fc Haehnelt I (|2000r ). It is possible that the 
correlation could be produced in this case if the opti- 
cal depth of the corona r increases with the mass of the 
black hole A/bh- This is equivalen t to assuming that £ 
increases with Mbh in the model of lHaardt et al.l ()1997f ) 
for an dominated corona. Since C oc Mbh, the previ- 
ous assumption is valid if we assume that the size of the 
corona is kept approximately constant as A/bh increases 
(see equation (HI ))- Under this assumption, both steady 
coron a models (|Haardt et al.l 119971 : IMerloni fc Fabianl 
l2001f) analyzed here will reproduce the F — Lx relation. 

The possible evolution of the slope and offset of the 
F — Lx correlation found in this study of RQ AGNs can 
be explained by an evolution of the properties of the hot 
corona. Specifically, the slope of the F — Lx correlation 
depends on the optical depth and compactness parame- 



ter of the corona in the model of Haardt et al. (1997) 
and the optical depth of the active regions in the thun- 
dercloud model by Merloni fc Fabian (2001). If AGNs 
within a certain redshift range contain hot coronae of 
similar properties we expect them to show a significant 
correlation between F and Lx in this redshift range. One 
explanation of the possible evolution of the slope and off- 
set of the F— Lx correlation is that the mean properties of 
the hot coronae of AGN evolve resulting in a detectable 
change in the slope and offset of the F — Lx correlation 
between AGNs at z ~ 2.2 and z ~ 0.7. One possi- 
ble explanation for the slight decrease in the strength of 
the F and Lx correlation in the second redshift bin is 
that this redshift interval is comprised of AGNs having 
a large range of coronal properties leading to a weaker 
correlation between F and Lx- We note that the possible 
decrease in the significance of the F and Lx correlation 
in the second redshift bin needs to be confirmed with a 
larger sample of RQ AGN. 

6. CONCLUSIONS 

In this paper we have selected a sample of radio-quiet 
AGNs (173) from the CDF surveys with moderate-to- 
high S/N, and have found strong evidence of a correla- 
tion between the X-ray spectral parameters F and Lx- 
We found that the slope and offset of a linear fit to 
the F — Lx correlation possibly evolves for sources with 
z > 0.1. Analyzing this relation in three different red- 
shift bins that contain a similar number of sources (~50) 
we conclude that this correlation is highly significant in 
two redshift bins, 0.3 S z < 0.96, and 1-5 < z < 3.3 and 
slightly less significant in the redshift bin 0.96 ^ z < 1.5. 
We note that the possible weakness of this correlation 
for sources with 0.96 < z < 1.5 appears to be driven by 
the absorbed sources in this redshift range. The F — Lx 
correlations in 0.3 < z < 0.96 and 1-5 < 2 < 3.3 are sig- 
nificant at the >99.9% confidence level for fits performed 
in the 0.5-8 keV observed- frame and at the >99.5% con- 
fidence level in the 2-10 keV rest-frame. 

The fact that this correlation is also present when we 
estimate the luminosities in the 2-10 keV rest-frame, and 
also holds for sources with low column densities, suggests 
that this correlation is not artificially driven by any un- 
modeled complexity in the intrinsic absorption (iVn). We 
performed several tests to investigate whether the F — Lx 
correlation found in this study is produced by a change 
with luminosity of the Compton-reflection component. 
We found that the strengths and slopes of the F — Lx 
correlation are similar for fits performed in the 0.5-8 keV 
observed-frame and 2-10 keV rest-frame bands. Our 
analysis indicates that the strengths and slopes would be 
significantly different if the correlation was driven by a 
Compton-reflection component. The difference between 
the observed weighted mean values of F obtained from 
fits performed in the 0.5-8 keV observed- frame and 2- 
10 keV rest-frame bands is less than 0.03. Our simula- 
tions indicate that if an un-modeled Compton-reflection 
component was producing the observed correlation a dif- 
ference of F of about 0.2 would be expected. We conclude 
that a Compton-reflection component is unlikely driving 
the F — Lx correlation found in this study. 

This correlation applies to sources with two different 
luminosity populations; one with log L2-10 ~ 43.1 ± 0.5 
(0-3 < z< 0.96), and the other with log L2-10 ~ 44.1 ± 
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0.4 (1.5 < z < 3.3) (see Figure [9]), indicating different 
populations of AGNs. The F — Lx relation results in 
a softening of the X-ray spectra as the luminosity of the 
AGNs increases. 

The F — Lx correlation found in the redshift range 
of 1.5 < z < 3.3 is of special interest because it con- 
firms a previous independent study of RQQ at z > 1.5 
pai et al.ll200l . Co mbining data f r omlDai et al.l (I2004D 
and other surveys (IVignali et al.l 119991: iGeorge et al.l 
120001: iReev es fc T urner I l2000t iPage et al.l |2003[ ). cited 
in iDai eTlil. (2004i) , we find that the F — Lx correla- 
tion becomes even more significant in the luminosity 
range of £2-10 ~ lO''^ - 8 x lO'^^erg s~^ with a Spear- 
man correlation coefficient of rg ~ 0.6 significant at 
the >99.9% confidence level (null hypothesis probabil- 
ity ~ 1.6 X 10"^). 

We presented two steady-cor ona models (jHaardt et al.l 
[19971 iMerloni fc Fabian llBoih that can explain both the 
F — Lx correlation found in this work and the saturation 
observed in the L2-10 vs. F relation using the surveys 
analyzed in §5.2.7. Based on these models, we proposed 
two different interpretations to explain the F — Lx corre- 
lation and its possible evolution with z. The first inter- 
pretation posits that this relation is driven by changes in 
the Eddington ratio (e = Lboi/-^odd) for a population of 
AGNs of similar mass. The second interpretation posits 
that the F — Xx relation is driven by changes in the mass 
of the AGNs. The present analysis does not allow us to 
infer which of these two scenarios is primarily responsi- 
ble for driving the F — ix correlation; however, future 
measurements of the black-hole masses for several AGNs 
in our sample will allow us to resolve this issue. 

To explain the detected possible evolution of the slope 
and offset of the linear fit to the F — Lx correlation we 
have proposed a simple model that posits that the mean 
properties of the hot coronae of AGN at z ^2.2 differ 
significantly from those of AGN at z ~ 0.7. This model 
also assumes that within each redshift bin the optical 
depths of the hot coronae of the AGNs are similar. 

We note that the detected change of the F — Lx corre- 
lation found in our study applies to RQ AGNs detected 
in the CDFs, which are two representative and normal 
fields. Further X-ray spectral studies of deep Chandra 
fields will test if the detected change of the F — Lx corre- 
lation also applies to wider fields of view. Expanding the 
sample will also test the correlation in narrower redshift 
bands and thus better constrain the epoch at which pos- 
sible changes in the average emission properties of AGNs 
occurred. 
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